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Abstract 

We give predictions for diffractive heavy flavour production at the Tevatron and 
the LHC in leading-order approximation. In the framework of these studies we use 
three different models for the partonic structure of the Pomeron recently proposed 
by Stirling and Kunszt. These Pomeron models are, despite being fitted to the 
same diffractive deep inelastic HERA data, very different in their parton content 
and taken together provide a powerful tool to probe the structure of the Pomeron. 
All models satisfy GLAP evolution and show a significant dependence. We give 
numerical predictions for single as well as double diffractive cross sections assuming 
a Donnachie-Landshoff-type Pomeron flux factor. 



1 Introduction 



Recent observations at HERA |T], 0] have provoked a renewed interest in an already long 
standing concept, the concept of the Pomeron in high energy scattering. This idea was first 
introduced by Pomeranchuk in 1958 to explain the fact that in the high energy limit 
the measured hadronic cross sections remain (approximately) constant. This asymptotic 
behaviour was explained in terms of a Regge trajectory a{\t\) with the Pomeron trajectory 
having the intercept a(0) ~ 1. 

Present and future generations of hadron-hadron and lepton-hadron colliders may be 
able to actually observe the Pomeron in diffractive events. We shall characterise an event 
as single diffractive if one of the colliding hadrons emits a Pomeron P that scatters off 
the other hadron. A single diffractive reaction would therefore correspond to 

hadron i + hadron j — > hadron i + X, (1) 

in which the intact hadron i emits a Pomeron (hadron i —>■ hadron i + W) which interacts 
with hadron j (hadron j + W —>■ X) leading to the reaction of Eq. (|l|). Hadron i is 
detected in the final state with a large longitudinal momentum fraction. Hard diffractive 
events, events with a large momentum transfer Q (1/Q -C 1 fm), are also characterised 
by a second feature: the absence of hadronic energy in certain angular regions of the final 
state phase space. These rapidity gaps, reported by the ZEUS and HI collaborations at 
HERA, are suggestive of a colour neutral object, the Pomeron, emitted from the incoming 
proton. Events fulfilling the conditions of large rapidity gaps and a (highly excited) hadron 
remnant in the final state are called single diffractive to distinguish them from those in 
which both colliding hadrons remain intact as they each emit a Pomeron, the so-called 
double diffractive events. The HI collaboration reported a fraction of about 6% of deep 
inelastic scattering events being single diffractive @]. The HERA data suggest that F 
has a hard structure. This observation has also been given support in an earlier UA8 
experiment. A partonic structure of the Pomeron was first proposed by Ingelman and 
Schlein in 1985 |^ to calculate high pr^iet production in single diffractive dissociation. 
The discovery of diffractive dijets (high p^-jets) with px > 8 GeV at the CERN pp-colhder 
in 1988 by the UA8 collaboration |^ fitted exactly with their prediction. 

In this paper we give predictions for single and double diffractive heavy flavour pro- 
duction at the Tevatron pp-collider {y/s = 1.8 TeV) and the LHC pp-collider, with a 
centre-of-mass energy ^/s = 10 — 14 TeV. By heavy flavour we mean QQ production 
with Q = c (charm), b (bottom) or t (top). Even though several tests of diffractive W, 
Z or Higgs production 0, ||, |, ^ have already been proposed, we shall focus especially 
on QQ states as their cross sections turn out to be typically larger and this makes their 
observation, especially in the diffractive context, easier in principle. We shall employ 
different aspects all deduced from collider physics phenomenology to create a solid model 
for diffractive scattering via Pomeron exchange. 
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The emission of a Pomeron from a fast-moving hadron will be treated in the framework 



of the Donnachie-Landshoff model of a reggeised, non-perturbative Pomeron. In the 
literature this type of Pomeron is called soft to discern it from the perturbative {hard) 
Pomeron. The number of Pomerons emitted by a hadron is described by the flux factor 
fv/i{x^i \ A) where i stands for any hadron. The longitudinal momentum fraction of the 
Pomeron and the t-channel invariant momentum transfer to the Pomeron fully determine 
the emission factor /p/i(a;^, |t|). As will be discussed in Section 2, a Regge factor is 
imbedded in the flux factor instead of a Pomeron propagator. 

The success of this prediction has been exploited to expand the ideas of the parton 
model to the Pomeron. The Pomeron, being a colour neutral object carrying the quantum 



numbers of the vacuum, should be made of qq pairs and/or gluons. Factorisation [|T2 
defines the diffractive structure function F.p{x,Q'^,x^, \t\) as a product of the structure 
function of the Pomeron F2{x, Q^) and the flux factor, i.e. 

F^{x,Q\x^M) = F^{x,Q') ® Mx^M). (2) 

Measuring the diffractive structure function F2 and assuming a particular flux (emission) 
factorS /F/i(a;^, |t|) gives information on the Pomeron structure function F2{x,Q'^) and 
hence on the parton content of the Pomeron. Because factorisation is assumed to give 
a universal picture of the Pomeron structure, the HERA data on F2{x, Q^) can be used 
to give predictions for hard diffractive scattering at hadron-hadron colliders. Comparing 
these predictions with measurements allows the factorisation hypothesis to be tested. 
This leads back to the fact that F^ has recently been measured at HERA |jl3| and 



thus provides new information on the partonic structure of the Pomeron. The data still 
suffer from large errors but nevertheless allow a fit of possible parton distributions inside 
the Pomeron. This has been done for example in Refs. p, ^ . In particular, in Ref . 



three Pomeron models were proposed to fit recent HI data [|1^. At this stage it has not 
been possible to decide experimentally whether the Pomeron is more gluonic or quark-like 
- deep inelastic scattering (DIS) only gives information on the quark content. Taking this 
deficiency into account, the three models range from pure glue up to pure quark content at 
the starting scale Qo = 2 GeV. While the quark content in each of the three models turns 
out to be rather similar, the gluon contents differ from soft to hard. No momentum sum 
rule for the distributions is imposed and an overall normalisation factor is incorporated 



into the flux factor, as described in Ref. |15|. A more detailed discussion can be found in 
Section 3. 

The aims of this paper are as follows. We want to show how the combination of Regge 
theory and factorisation can be applied to calculate diffractive cross sections, especially 
in the light of recent experimental F2 data. The theory of this is reviewed in Section 2. 
So far the theory has been mostly applied to single diffractive events. The reason is that 



^The flux factor can also be extracted from experiment. 
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double diffractive cross sections at present and future colliders are generally too small to 
play a pivotal or even measurable role. But especially after experiences with top quark 
production and recent measurements of charm and bottom quark total cross sections 
at the Tevatron [|I^ and elsewhere, it is worth looking for heavy flavour production in 
both single as well as double diffractive cross sections. Future measurements at high 
energy hadron-hadron colliders together with our predictions could be a powerful step 
towards a better understanding of the Pomeron structure, in particular to reveal a gluonic, 
quark-like or mixed Pomeron. We shall also consider the next-generation collider, the 
LHC, for this purpose. With a centre-of-mass energy ^/s = 10 TeV, we shall show that 
double diffraction will be no longer an unmeasurable effect but a paradigm for visualising 
Pomeron-Pomeron scattering. 

We note that the concept of a universal factorisation of Pomeron structure, in partic- 
ular for hadron-hadron collisions, has recently been questioned, see for example Ref . jl^ . 
Our predictions and further experiments at the Tevatron and the LHC might help to 
clarify this question. 

All numerical predictions for the Tevatron and the LHC are presented in Section 4. 
The results are critically discussed and the different models are carefully discerned. Both 
absolute values as well as single and double diffractive ratios are shown. Section 5 sum- 
marises the results, underlines the most important features and discusses open problems. 



2 Kinematics and Cross Sections 

In this paper we concentrate on the single a^^ and double diffractive a^^ cross sections 
for the production of heavy flavour QQ states. What do we mean by diffractive? First let 
us consider the case of single diffraction, where a Pomeron with momentum kf is emitted 
by one of the colliding hadrons as shown in Fig. la. The definition of single diffractive 
means: one hadron, emitting the Pomeron, is being detected (at least in principle) in 
the final state. The other hadron scatters off the emitted Pomeron. A typical single 
diffractive reaction for our purpose at the LHC would read: p + p — > p + QQ + X. The 
Pomeron kinematical variable xf is defined as: xf = sj/ Sij , where is the centre-of- 
mass energy in the Pomeron-hadron j system and y^Si] = \fs the centre-of-mass energy 
in the hadron z-hadron j system. Instead of the gluon- or (anti) quark-distributions of 
hadron i {fg,q/i{xi,Q'^)), the parton distributions of the Pomeron {fg^q/w{xi/xf,Q'^) with 
Xi < xf) emitted by hadron i define the single diffractive heavy quark cross sections 
assuming a partonic picture of the Pomeron. We have therefore to substitute 

/f dx^ 
d\t^\ j fg,g/w{Xi/xl Q^) /p/i (xf , \ti\) . (3) 

\U\ 

1 
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The Pomeron parton distribution fg,q/w{xf,Q^) is weighted by a flux factor fw/i{xf, \ti\) 
which gives a measure of the number of Pomerons emitted by hadron i. It is therefore 
a function of the energy loss of the proton tj = {Pi — P/)^ due to Pomeron emission as 
well as of the momentum fraction xf carried by the Pomeron. Its form can be modelled 



by 111 



/i>Mxf,|t.|) = ^[i^i(|t.|)f(xf) (4) 

The coupling of the Pomeron to the proton is assumed to be 3/5-Fi(|ti|) with /3 = 1.8 GeV~^. 
The factor 3 corresponds to the number of valence quarks inside the hadron. This is dis- 
cussed and supported in |n|] . In fact it should couple in a pointlike form to single quarks. 
From this point of view the effective Pomeron structure function should contain quark- 
antiquark contributions. An analogy with the photon structure function emerges from 
this picture. The elastic form factor of the proton is experimentally determined 

and well parametrised by |Tl 

4M,^ + 2.8|t.| / \U\ y 

where Mp is the proton mass. Instead of propagators for the Pomeron, one inserts the 
Regge factor (^xf^ '''' with the Regge trajectory = 1 + e — a'\ti\ with e and a' 



both deduced experimentally. Latest experiments at the Tevatron [|T8l showed a slight -^i 
dependence of e as long as a' is fixed. The CDF collaboration studied single diffractive 
events p + p ^ p + X and obtained from a data fit: 

e = 0.121 ±0.011 at = 546 GeV, 
e = 0.103 ±0.017 at = 1800 GeV. 

They conclude that large screening effects have to be introduced to save the traditional 
Pomeron model. An analysis at ^/s = 20, 546 and 1800 GeV showed a significant 
^0.030 dependence for the single diffractive cross section. 

In practice, e, introduced as an effective power, measures the deviation of the differ- 
ential cross section da^^ / dM"^ from a pure 1/M^ dependence. On the other hand, taking 
only the triple-Pomeron diagram as contribution to diffractive dissociation into account 
and neglecting screening effects, 1 ± e can be interpreted as the intercept of a supercritical 



Pomeron 19 



Also deduced experimentally by CDF was the sensitivity to a' . A change in a' by 
5a' = ±0.1 GeV~^ results in a change in a^^ of only ±0.1% and e changes at the same 
time by Se = ±0.011. Neglecting any s dependence, we can fix e and a' to the standard 
values |ni e = 0.085 and a' = 0.25 GeV'^. 



According to [11 all these relations should only be valid if we limit the Pomeron 



momentum fraction to xf < 0.1. Thus only a maximum energy loss of 10% for each 
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hadron is allowed. But assuming the final state proton is not detected, there will be no 
energy cut for t, although the flux factor decreases very fast with \t\. 

The flux factor plays a pivotal role because it directly yields the energy dependent 
emission rate of the Pomerons and is therefore extensively discussed in the literature. 
Bruni and Ingelman calculated the diffractive W production rates at the Tevatron using 
a constant Pomeron-proton cross section of 2.3 mb, obtained from Regge analysis but 
with a non-Regge parametrisation for the flux factor ||^ 

Mxl \U\) = ^(6.38e"«l*'l + 0.424e-=^l*»l)i-, (6) 

z.o 

fitted from single diffractive data. Throughout this work, however, we shall use the 
Donnachie-Landshoff Pomeron flux factori of Eq. . 

We shall also consider double diffractive events. The procedure is analogous to the 
single diffractive scenario and is sketched in Fig. lb. Now both hadrons undergo Pomeron 
emission and therefore all hadronic parton distributions are replaced by the parton dis- 
tributions of the Pomeron via the rule given in Eq. (|^). However, the basic formulation 
of diffractive scattering was intended for reactions where one hadron scatters diffractively 
and one hadron is highly excited. The crucial point is that the single diffractive cross sec- 
tion of a hadron-hadron collision is assumed to factorise into the total Pomeron-hadron 
cross section and the Pomeron flux factor 0. The single diffractive event may then be 
written as 

^^P^l^ I = h/iiXi ,\ti\)(8)a ( F+ p{p) ^QQ + X). (7) 

Eq. is the analogue of Eq. (H). Factorisation has also been applied to double diffraction, 
for example predictions for Higgs |jlO[ or W boson production p, |^ at the Tevatron and/or 



the LHC Both colliding hadrons can in principle be detected in the final state. Again, a 
typical reaction at the LHC would look like: p + p — > p + p + QQ + X. Double diffractive 
events thus are characterised by two quasi-elastic protons with rapidity gaps between 
them and the central heavy flavour products. 

Using the substitution given in Eq. (^) we can proceed in complete analogy to the 
single diffractive case and finally write down the expressions for the total, single and 
double diffractive cross sections for QQ production 



1 1 



a'^{s,mQ) = j j dxidx2 fq/i{xi,Q'^)fq/2{x2,Q'^) X aqq{s,mQ,as{Q'^)) 



X\=T X2 = — 

^1 



^Goulianos has recently argued |20 that only a normalisation of the flux factor to one Pomeron per 
nucleon might yield physical results. 
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+ fg/iixi,Q^)fg/2ix2,Q^) X agg{s,mQ,asiQ^)) 



a 



SD 



dxi / dx 



d\ti\® U/i{x\,\ti\) 



a 



DD 



{s,mQ) 



X [fq/¥iXl/x^,Q'^)fq/2iX2,Q^) X agq{s , TUq , OsiQ^)) 
+ fg/Tp{xi/x^, Q^)fg/2{.X2, Q^) X agg{s, TUq, OslQ^)) 

+ (1^2), 

dx^ 



(9) 



dxi / dx 



x^ 



(J/X'2 



d\ti\ / d\t2\®hii{xlAti\)h/2{xl\t2\) 



/g/F(xi/a;f,(5^)/g/ff(x2/x2,Q^) x 0-5^(5, mg, as(Q^)) 
+ fgiw{.Xi/x\,Q'^)fg/^{x2/xl,Q'^) X agg{s,mQ,as{Q'^)) 



(10) 



Of course the threshold condition s = X1X2S > ts = 4mQ for the production of a QQ pair 
of mass 2mQ must be fulfilled and the cut-off in the Pomeron spectrum Xj < xf < 0.1 
taken into account. With s we define the Mandelstam variable of the two contributing 
subprocesses qq QQ and gg — » QQ and r = 4mQ / s regulates the threshold condition 
automatically as function of the integration boundarieg^l. We restrict ourselves to the 
leading-order contributions. The Feynman diagrams of this Born approximation are 
shown in Fig. 2. A discussion of next-to-leading-order non-diffractive heavy flavour 
production can be found in [pl|] . 

We only cite the expressions for the subprocess cross sections in leading-order approx- 
imation. With p = ^toq I s they read 



gg 



(J. 
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's,mQ,as{Q^)) 



[s,mQ,as{Q'^)) 



1 7ia\Q^ 
48 s 



(112 + 124p)Jl-p 



+ (16 + 16p + p^) In 



1 + 



levra^ 
27 I~ 



(11) 
(12) 



where the renormalisation scale is taken to be the subprocess collision energy, Q"^ = s = 
X1X2S. For an analytic derivation of Eqs. (11) and (12) using standard techniques we refer 
the reader to l2g,|2|,g. 

One might however argue that for the gg —>■ QQ subprocess the t- and u-channel con- 
tributions should dominate [^], in which case a choice of = | (("^g — i) + {mq — u)^ 

= |s might be more reasonable. In fact the parton distributions are affected by the 
choice of Q^ in the framework of GLAP evolution. But as in [E3] we could not find a 



For the sake of simplicity we did not explicitly write down the integration boundaries for a and 
a^^ . They can easily be deduced from the two conditions: (1) X1X2 > t and (2) Xi < xf < 0.1. 
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significant sensitivity of our results to tliis clioice compared to a general Q"^ = s for botli 
subprocesses. 

A still unsolved problem is the structure of the Pomeron. Nearly all possible scenarios 
of parton distributions proposed have favoured a more or less gluonic content. In the 
following section we shall give a very cursory presentation of the most popular Pomeron 
models at this stage. A pivotal role will be played by three Pomeron models, recently 
presented and applied to diffractive W production at the Tevatron and the LHC by 
Kunszt and Stirling They will be the input distributions for the diffractive heavy 
flavour production, which we shall deal with in Section 4. 



3 Models of the Pomeron 

The Pomeron has been postulated to describe hard diffractive collider phenomenology. 
In the framework of Regge theory a colourless object, carrying the quantum numbers of 
the vacuum, is able to explain the new observations. A pure gluonic composite seemed to 
be the simplest explanation on the basis of partonic contents and proposed features ||25 |. 



The Pomeron was therefore assumed to behave essentially as a hadron and Ingelman and 
Schlein introduced the concept of a Pomeron structure function . 

Further collider experiments, by for example the UA8 collaboration at CERN ||^, gave 
evidence for a hard parton distribution inside the Pomeron but still could not distinguish 
between a gluon or quark dominated Pomeron. 

Another striking problem, not yet fully resolved, was the character of the Pomeron: 
should it be treated by perturbative or non-perturbative means? Landshoff argued 



that experimental results are best described by a Pomeron having propagators that are 
non-perturbative. They must not have a pole at /c^ = 0. He summarises the phenomeno- 
logical features of the so-called soft Pomeron as follows: 

• it is rather like a C = +1 isoscalar photon, 

• it couples to single quarks with a strength /? = 1.8 GeV^\ 

• it is a simple Regge pole with a Regge trajectory a;(|t|) = 1 + e — a'\t\ taking into 
account an effective power e that decreases with increasing ^/s, but very slowly such 
that e ~ 0.08 at typical hadron collider energies. 

Studies at the pp-collider at CERN favoured a more hard quark- or gluon-like 
structure function for the Pomeron than a soft onea. Even though the Pomeron exchange 
is a pure gluon exchange, this does not imply that the Pomeron structure function is 
entirely, or even predominantly, gluonic. 



''The conclusion of UA8 was that ^ 57% of aU events should be hard \Qx{l — x)] and ^ 13% were 
tagged to be soft [6(1 - xf]. The remaining 30% showed a superhard [5(1 — x)] characteristic. 
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This was the starting point to give predictions for diffractive events using different 
Pomeron models, for example for diffractive W and Z boson production or for diffrac- 
tive Higgs production . 

Whatever the parton distributions looked like, in the original investigations they were 
assumed to show no evolution with the scale at which the Pomeron is probed. They 
were functions only of the momentum fractions of the partons inside the Pomeron. This 
situation has changed, since there are now data for the diffractive structure functions 
measured at HERA |^, Q at different values of Q^. Employing Eq. and assuming 
factorisation directly yields data for the Pomeron structure function F2^(x, Q^), assuming 
a certain structure for the flux factor. This has been discussed in detail in the Introduction. 

We shall focus on the latest set of three different Pomeron models, based on recent 
HERA data [|l^, provided by Kunszt and Stirling The main features of these models 
are: 



• all three models are qualitatively very different concerning the partonic contents at 
the starting scale Qo = 2 GeV, 

• all three models give satisfactory agreement with the HI and ZEUS data, 

• all models undergo leading-order GLAP evolution as for the usual parton distribu- 
tions. 

In qualitative terms, the three models can be characterised as: 

Model 1: At Qo the Pomeron is entirely composed of quarks. Gluons are dynamically 
generated via GLAP evolution; 

Model 2: A mix of quarks and gluons at starting scale Qq; 

Model 3: A predominantly hard gluonic content at starting scale, the gluons inside the 
Pomeron carry large fractional momenta. 

An SU(3) flavour symmetry is assumed for all three models, such that a momentum 
sum rule for Qq = 2 GeV can be imposed 

1 

J dxx (Qfq/Tpix,Ql) + /g/F(a;,Qo)) = 1, (13) 



to deduce the overall normalisation factor for the data fit. Because the quark contents 
inside the three models are more or less comparable since they are constrained by the 
HERA data, we show in Fig. 3 solely the Q evolution of the gluon distributions xfg/^^x, Q"^) 
for the three Pomeron models. 



8 



We shall see later that it is the gg contributions which predominantly govern the 
behaviour of the (diffractive) QQ cross sections for the Tevatron and the LHC in the 
given energy regime. These cross sections can therefore in principle provide important 
information on the gluon content of the Pomeron, complementary to the constraints on 
the quark content from DIS structure functions. For a broader discussion of the partonic 
Pomeron contents and their fit to the HI data we refer the reader to [^. 

It is very important to have distributions that undergo GLAP evolution, because, 
ranging from charm to top quark masses, different scales are probed and, as can be 
seen in Fig. 3, the distributions evolve significantly with Q^. This is different from most 
of the former static distributions proposed and enables us to give predictions for such a 
large mass range covered by the heavy flavours. 

The main features of the three Pomeron models concerning xfg/^{x, Q"^) can easily be 
deduced from Fig. 3. The left side shows surface plots of the gluon distributions in the 
regime 4 GeV < Q < 50 GeV and 0.1 < x < 1.0. The corresponding contour plots in the 
x-Q plane are presented on the right side. 

As can be seen in Fig. 3, Model 1 gains only a small gluon distribution via GLAP 
evolution at Q = 4 GeV. Recall that there are no gluons at all at the starting scale 
Qo = 2 GeV. The gluon distribution is almost flat in Q for moderate and high values of 
X. This fact can also be seen in the contour plot, where for Q > 10 GeV there is only a 
significant evolution in the very small x regime. No gluons with x > 0.6 are present in 
this regime. 

This is different to Model 2, which contains gluons with moderate or even higher 
fractional momenta. The gluon distribution shows a significant Q evolution not only for 
small but also for moderate x (this can be seen clearly in the contour plot), even though 
the distribution asymptotically approaches zero for higher Q and x > 0.3. 

Model 3 finally shows the most dynamical evolution of the gluon distributions not only 
in Q but also in x. Very hard gluons are present for high values of x. This is due to the 
basic construction of Model 3 as discussed above. These gluons lose their high momenta 
at larger values of Q mainly due to quark-antiquark pair production. The depletion of 
these high fractional-momenta gluons with increasing Q again becomes very transparent 
in the corresponding contour plot. For a further discussion we again refer the reader to 



4 Diffraction at the Tevatron and the LHC 

For the GLAP evolution of the parton distributions inside the (anti)proton we use the 



MRS (A') set presented in One test of the reliability of this set is the comparison of 
the calculated total inclusive cross section for ti production at =1.8 TeV with the 
latest experimental measurement The CDF collaboration measures a total inclusive 
cross section of crj^p ip + p ^ ti + X) = 7.6l|o pb for an experimentally determined top 
mass mt = ( 176 ± 8(stat.) ± 10(syst.) ) GeV, while the MRS(A') set of partons predicts 
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a'^{p + p ^ tt + X) = 5.81 pb for rrit = 176 GeV. Measurements of inclusive cc and bb 
cross sections at the Tevatron collider yield ||16| : 

^Jxp.(P + P^cc + X) = 115.2^24-3 nb, 
aljp + p^bb + X) = 30.1™ nb 

for rric =1.5 GeV and = 4.5 GeV compared to the leading-order theoretical predictions 
of a'^{p+p cc+X) = 98.7 nb and cr^{p+p bb+X) = 44.2 nb for the same and rrif,. 
Even though we are working only in leading-order approximation, the agreement with 
the experimental results is encouraging. Predictions for the ratios a^^/a'^ and a^^/a'^, 
however, should be relatively insensitive to higher-order perturbative corrections. 



4.1 Probing the Pomeron at the Tevatron 

We calculate numerically the total inclusive, the single diffractive and the double diffrac- 
tive cross sections for heavy flavour pair-production using the MRS (A') set of partons 
and the three Pomeron models introduced in p and discussed in Section 3. Our results 
are shown in Fig. 4. Supplementary we show the pure gluon fusion contribution as dashed 
lines. 

For the Tevatron collider the qq process becomes dominant for tuq > 50 GeV. This qq 
dominance is obviously visible in the total cross section but becomes even more striking 
for the single diffractive one. Due to this process, the single diffractive ratio R^^ defined 
as even increases for mq > 70 GeV before this process runs out of 

centre-of-mass energy due to the cut-off in the Pomeron spectrum {xf < 0.1) and 
steeply falls to zero. Note that this restriction on a/^ unfortunately takes place before 
the top quark production domain is reached. The single diffractive ratio R^^ for all three 
Pomeron models reaches its local maximum at slightly lower masses than the top quark 
mass. This gives little hope to observe diffractive top quark events at the Tevatron but 
is promising for the LHC as will be discussed in Section 4.2. 

At Tevatron energy the threshold for single diffractive events is reached when mq = 
270 GeV. But single diffractive events might be observed only up to mq ^ 150 GeV with 

Model 3 : a^^^ = 1.05 pb ~ 4.05% of , 
Model 1 : a^^f^ = 0.22 pb - 0.86% of 

for niq = 150 GeV. This is qualitatively what we expect. Model 3 with the hard gluon 
content gains more and more qq pairs via GLAP evolution which dominantly contribute 
to the single diffractive cross section for higher masses. Model 2 with the quark-gluon 
mixture dominates for small masses, but at mq ~ 10 GeV Model 3 takes over. The 
production rate of the important qq pairs for higher masses increases more rapidly in 
Model 3. Model 1, however, without any gluons at the initial Qq = 2 GeV scale, marks 
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the lower limit for single and double diffractive scattering over the whole mass range. The 
numerical values for the single and double diffractive production of charm, bottom and 
top are given in Table 1. 

Let us now focus on Models 2 and 3: the ratios for single diffractive scattering range 
between ~ 22% for charm and ~ 15% for bottom quark production and reach their 
locally lowest values of ~ 2% for uiq = 70 GeV, before the ratio even increases due to 
the qq contribution in the mass region 70 GeV < mg < 150 GeV locally peaking at 
ruQ = 150 GeV with ratio R^^ = 5.5%(4.2%) for Model 3 (2). Qualitatively Model 1 
shows the same behaviour. But with only quarks in the starting distribution, the creation 
of gluons and qq pairs via GLAP cvohition proceeds only slowly. This can be seen in 
Fig. 3. The difference is expressed in the large gap between Model 1 on the one hand and 
Models 2 and 3 on the other hand. So the conclusion is that at least for single diffraction 
a gluon rich Pomeron as input should be more easily detected. Model 1 and Model 3 for 
example differ by about one order of magnitude. This turns out to be a crucial difference 
with such absolutely small diffractive cross sections. 

Double diffractive scattering seems to favour a balanced mixture of quarks and gluons 
in the starting distribution and during GLAP evolution, as provided in Model 2. This 
can again be observed in Fig. 4. 

Notice that the former dominant hard gluon Model 3 quantitatively shows about the 
same behaviour as Model 1. Why the behaviour of Model 3 in double diffraction is 
different from that in single diffraction can be understood by a further analysis of the 
kinematics among the partons inside the Pomeron and their distributions at different 
mass scales. Models 2 and 3 show a crossing in the case of single diffractive scattering 
as can be observed in Fig. 4. Its existence can be immediately explained in terms of 
the gluon distributions which are shown in Fig. 3. As we have already pointed out, the 
gluon distributions govern the behaviour of the cross sections at the Tevatron, especially 
for small quark masses. For small x and Q, the gluon distribution of Model 2 is slightly 
bigger than that of Model 3. For higher x and/or Q the situation is reversed. 

In the case of double diffraction, an analysis of the average fractional gluon momentum 
{x) = (xi/xf) yields (x) ~ 0.11 for mg = 2 GeV and (x) ~ 0.24 for tuq = 50 GeV for 
all three models. In this regime the gluon distribution of Model 2 again exceeds that of 
Model 3. So, no crossing can be observed and Model 2 dominates throughout. In fact 
the crossing would take place at mg ~ 80 GeV, shortly before the threshold for double 
diffraction with mg ~ 90 GcV is reached. 

The relatively broad gap between Model 2 and Model 3 appearing in Fig. 4 can again 
be explained by the same straightforward analysis of the average gluon momentum. Again 
a comparison with the corresponding gluon distributions of Model 2 and Model 3 in Fig. 3 
shows the absolute difference of these distributions in the regime 0.1 < (x) < 0.2. 

While diffractive charm and bottom quark production might be observable at the 
Tevatron, there is no hope for diffractive tt pair production being visible at a total centre- 
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of-mass energy of i/s = 1.8 TeV. Either the top mass exceeds the kinematic threshold 
(double diffraction) or the effect of the threshold is already strongly influencing the process 
by a steep decrease in the diffractive cross section (single diffraction near top mass). 
However at the LHC such heavy flavour threshold suppression is less severe, as we shall 
see in the following Section. 

4.2 Predictions for the LHC 

The LHC will provide a rich field of study for diffractive events, even in the top mass 
regime. With a centre-of-mass energy of at least y/s = 10 TeV |^ the double diffractive 
threshold is reached for ttiq = 500 GeV and the single diffractive one lies at rriQ = 
1500 GeV. This upper bound is very promising in particular for diffractive top quark 
production at the LHC. 

At a proton-proton collider, the dominant process for QQ production is of course 
gluon-gluon fusion since antiquarks only appear as sea quarks inside the proton. For a 
top mass of 176 GeV the pure gg contribution is about 91% of the total cross section in 
our calculation. This also holds for the single and double diffractive case. 

The numerical results are shown in Fig. 5 and the numerical values are again listed 
in Table 1. The single diffractive ratios R^^ for the three Pomeron models are between 
20 — 40% for charm and 10 — 40% for bottom quarks. The maximal and minimal single 
diffractive top quark rates are 



with rrit = 176 GeV. Even though the single diffractive ratios for top production are 
comparable to the Tevatron rates, the absolute single diffractive cross sections are crucially 
enhanced. We find an enhancement of about a factor 100 at the LHC. For example Model 2 
and Model 3 provide a single diffractive cross section of approximately 5 pb. This is about 
the total cross section for top production at the Tevatron. So, single diffractive top quark 
events should be readily detected at the LHC. 

The predictions for double diffractive scattering, however, are still not very promising. 
The maximal double diffractive cross section (for Model 2) is 1.86- 10^^ pb for top quarks. 
The qualitative behaviour of the three different Pomeron models is the same. Because the 
centre-of-mass energy at the LHC is larger by a factor of six, one can conclude from the 
qualitatively similiar behaviour of all three models at the LHC that the gluon distributions 
in all three models do indeed become similiar at higher Q, as already seen in Fig. 3. 

All conclusions that were drawn for the Tevatron still hold for the LHC, except that 
because of the higher centre-of-mass energy there are no kinematical artifacts in the 
considered flavour mass regime. Even for the top quark mass, the single as well as the 



Model 3 : cr; 



.3D 



5.45 pb ~ 1.66% of a'^, 
1.13 pb^ 0.34% of 



max. 
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double diffractive cross section behave rather smoothly. The influence of the threshold 
does not seriously affect the cross sections in this case. 

Considering the ratios for single and double diffractive events in Fig. 5 again shows 
a qualitatively comparable picture to the Tevatron. Again an analysis of the average 
fractional gluon momentum {x) will explain the differences in R^^ and R^^ . For R^^ , 
Models 2 and 3 are quantitatively equivalent, especially for the charm and bottom quarks. 
For a bottom quark mass of mi, = 4.5 GeV, we obtain as average fractional gluon mo- 
mentum (x) = 0.42 for Q ~ 2mj, = 9 GeV in the single diffractive case. A comparison 
with the gluon distributions in Fig. 3 shows that they are roughly equal for Model 2 and 
Model 3 in this region of x. This fact becomes even more transparent in the contour plots 
of Fig. 3. 

For the case of double diffraction the same analysis yields a lower average fractional 
momentum for the gluons, due to the energy-cut among both Pomeron emitting hadrons. 
Again for = 4.5 GeV we obtain (x) = 0.22. But in this region of x, the gluon 
distribution of Model 3 shows a local minimum, the hard gluons in this model give a 
rise of x/g/p(x, Q"^) only for x > 0.5 in the low-Q regime. The gluon distributions inside 
Model 2 also show a local minimum around x ~ 0.2, but its absolute value is higher than 
that for Model 3 in this region. This fact is responsible for the gap between Model 2 
and Model 3 as observed for R^^ in Fig. 5. For higher quark masses (higher values of 
Q'^) all three models become comparable concerning the gluon distributions, as already 
discussed. 

Model 2 and Model 3 as descriptions of the parton distributions of the Pomeron yield 
very promising single diffractive ratios, at least for charm and bottom quarks. Model 3 
predicts a single diffractive ratio of ~ 40% for cc and bb production. This is quantitatively 
comparable to the predictions of Model 2 as can be observed in Fig. 5 and numerically 
verified in Table 1. Thus about one third of the production of heavy flavours at the 
LHC including charm and bottom quarks should be single diffractive. But even Model 1, 
purely quark-like at starting scale Qo = 2 GeV, gives a single diffractive contribution of 
approximately 20% in this mass regime. 

5 Conclusions 

We have presented predictions for diffractive heavy flavour production at the Tevatron and 
the LHC. The input Pomeron distributions |Q] were different in contents from pure quark 
starting distributions up to hard gluonic contents at the GLAP starting scale Qo = 2 GeV. 
Our predictions can therefore be considered as upper and lower limits for the diffractive 
cross sections. 

We found that at the Tevatron single and double diffractive charm and bottom quark 
production is observable with a single diffractive ratio R^^ that is between 7% (charm) 
and 4% (bottom) for Model 1 as a lower limit and 26% (charm) and 13% (bottom) for 
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Model 2 as an upper bound. The double diffractive cross section for charm and bottom 
production lies above the total cross section for the top quark and might also be observable. 
The lower limit is governed by the pure quark-like Model 1 with 45 pb for charm 

and 18 pb for bottom quarks, whereas having Model 2 as the underlying Pomeron 

model gives a double diffractive cross section which is enhanced by a factor of four. 

The diffractive cross sections for charm and bottom quarks will be larger by three 
orders of magnitude at the LHC, and also single diffractive top quark production should 
be observable. The single diffractive tt cross sections range from 1 pb for Model 1 up to 
5 pb for the hard gluon Model 3, which is roughly the total top quark cross section at the 
Tevatron. All predictions are collected in Table 1 and shown in Figs. 4 and 5. 

The three Pomeron models presented in P] and used in our calculations are only a 
guide to possible Pomeron structures. Their determinations have been made possible 
by two strong assumptions: first, the Pomeron has a partonic structure and second, we 
assume factorisation such that recent measurements of the diffractive Pomeron structure 
function give a universal picture of the Pomeron. Only the latter makes predictions for 
the hadron-hadron collider meaningful. 

Our predictions can only be a guide as the calculations are illustrative only. We need 
a more detailed study with detector acceptance. This will follow as soon as there are 
first results on diffractive heavy flavour production at the Tevatron. But a study of these 
events is worthwhile, at least for diffractive charm and bottom production. Diffractive top 
quark events still will require a higher centre-of-mass energy before they can be studied. 
The LHC will be such a laboratory for diffractive scattering studies. We predict that 
about 30% of all charm and bottom production events should be single diffractive. 

It is not yet certain whether the theoretical framework of the Donnachie-Landshoff 
model can be extended to double or even multiple Pomeron exchange. Again, we give 
predictions for the case of an analogue treatment of double Pomeron exchange that will 
be interesting to test at future hadron-hadron colliders. Also the parameters used in 
this theoretical model have to be fitted to new high energy collider experiment J^. Finally 
future high energy experiments will hopefully more clearly reveal the structure of the 
Pomeron and answer some questions that were raised in this paper. 
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Table 1: The values for single and double diffr active cross sections, as well as their ra- 
tios to the total cross sections, are shown for average quark masses {rric) = 1.3 GeV, 
(mj) = 4.3 GeV and {rrit) = 176 GcV. The three different Pomeron models arc discussed 
in Section 3. We obtain numerical data for both the Tevatron and the LHC. A Q indicates 
that the threshold for this process was exceeded. 
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a) b) 

Figure 1: Kinematics of the single (a) and the double diffractive scattering (b) events, 
leading to QQ pair production. The subprocesses taken into account are either qq — > QQ 
or gg-^QQ . 
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Figure 2: The leading-order Feynman diagrams for the subprocesses qq — > QQ (a) and 
gg^QQ (b). 
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Figure 3: The gluon distributions of the three different Pomeron models are shown, 
starting with Model 1 at the top. For each model we show the surface plot for a given 
x-Q regime as well as the corresponding contour plot in the x-Q plane. 
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Figure 4: DifFractive Scattering at the Tevatron: The numerical resuhs of the total, 
single and double diffractive cross sections for the three Pomeron models are shown. 
The solid lines indicate the subprocess gg + qq ^ QQ, the dashed lines show solely the 
contribution from gluon fusion {gg — > QQ): (a) gives the absolute numbers for the cross 
sections (in pb), and (b) gives the ratios {R^^ = a^^/a'^ and R^^ = a^^/a'^). The 
mass regions of the charm (c), bottom (b) and top (t) quarks are indicated. The fixed 
centre-of-mass energy is 1.8 TeV. 
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Figure 5: Diffractive Scattering at the LHC: Same as Fig. 4 but now for a cen- 
tre-of-mass energy of 10 TeV. The sohd hnes indicate the subprocess gg + qq ^ QQ, the 
dashed hnes show solely the contribution from gluon fusion (gg — > QQ): (a) gives the 
absolute numbers for the cross sections (in pb), and (b) gives the ratios = cr^^/cr^ 
and R^^ = a^^/a^). 
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